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Abstract 

We construct the first explicit, smooth, horizonless black-hole microstate geometry 
whose moduli space is described by an arbitrary function of one variable and is thus infinite- 
dimensional. This is achieved by constructing the scalar Green function on a simple D6-D6 
background, and using this Green function to obtain the fully back-reacted solution for a 
supertube with varying charge density in this background. We show that this supertube 
can store parametrically more entropy than in flat space, confirming the entropy enhance- 
ment mechanism that was predicted using brane probes. We also show that all the local 
properties of the fully back-reacted solution can, in fact, be obtained using the DBI ac- 
tion of an appropriate brane probe. In particular, the supergravity and the DBI analysis 
yield identical functional bubble equations that govern the relative locations of the centers. 
This indicates that there is a non-renormalization theorem that protects these functional 
equations as one moves in moduli space. Our construction creates configurations that are 
beyond the scope of recent arguments that appear to put strong limits on the entropy that 
can be found in smooth supergravity solutions. 
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1 Introduction 



One of the key aspects to understanding black holes within string theory involves determining 
the families of supergravity solutions that have a given set of charges and angular momenta, and 
examining the relationship between the number of these solutions and the Bekenstein-Hawking 
entropy of a black hole with identical charges and angular momenta. 

The interest in this question is spurred by the desire to establish if, and in what form, the so- 
called fuzzball proposal (see [H El El HI |5] for reviews) applies to BPS black holes. This proposal 
requires that, within string theory, there should exist a huge number of horizonless configurations 
with unitary scattering, that have the same mass, charges and angular momenta as a black hole, 
and that start differing from each other at the location of the would-be black hole horizon. Fur- 
thermore, the number of such configurations should reproduce the Bekenstein-Hawking entropy 
of the black hole. If enough such configurations exist, then the AdS-CFT correspondence would 
strongly suggest that the classical black hole solution only gives a thermodynamic approximation 
of the physics that is correct at scales larger than the horizon scale, but does not correctly de- 
scribe the physics at the horizon, much as the thermodynamic description of a gas breaks down 
at scales smaller than the mean free path. 

The construction and investigation of black hole microstate configurations that has been going 
on over the past few years has yielded several interesting pieces of information. First, there is 
a very interesting geometric transition that takes the original black-hole or black-ring geometry 
over to a huge moduh space of smooth, horizonless bubbled geometries O [3 12] • Amongst such 
bubbled geometries there exist a very large number of smooth horizonless solutions that have the 
same charges and angular momenta as a BPS black hole with a macroscopically-large horizon 
area, both in five and in four dimensions |8l |9l [10]. Furthermore, these solutions can have a very 
long throat, and one can argue [HI [H] that they should be dual to CFT states that live in the 
same CFT "typical" sector as the states that give rise to the Bekenstein-Hawking entropy of the 
black hole. 

For practical, computational reasons it is easiest to study the microstate solutions that have a 
Gibbons-Hawking (GH) base space and thus descend to multi-center solutions in four dimensions 
jl2[ [13] . There are very large families of such solutions, however the semi-classical quantization of 
the moduli space shows that these solutions give an entropy that is parametrically smaller than 
the entropy of a black hole with the same charges |TTl jll] . This is to be expected, and indeed 
was anticipated in the earlier, very limited estimates of [8] : One cannot hope that counting only 
microstates that have a large amount of symmetry could reproduce the entropy of a given system. 
The simple bubbled solutions in five dimensions have the symmetries of the compactification torus 
as well as a tri-holomorphic U{1) isometry along the Gibbons-Hawking (GH) fiber. From such a 
class of solutions one cannot even reproduce the entropy of the two-charge D1-D5 system, which 
is known to come from smooth horizonless solutions [151 [13 [U] • The entropy of supergravity 
[/(l)-invariant solutions was also matched against the entropy of a weakly coupled graviton gas 
in AdS^ X S"^ in [14j. While this entropy is not restricted to [/(l)-invariant configurations, it is 
intrinsically perturbative, and cannot capture the supertubes we consider here; in the D1-D5-P 
frame our configurations give rise to a non-trivial topology, and in other frames have dipole 
charges corresponding to non-perturbative brane configurations. 

To some extent, the amazing part of f/(l)-invariant bubbling geometries and of four- 
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dimensional multi-center solutions is not the fact that there are so few microstates with a set 
of isometries (in the D1-D5 system there is only one such microstate), but that there are so 
many. The importance of the bubbled GH geometries is not so much in the large number of 
highly symmetric solutions but in the fact that they come from a geometric transition of branes 
in flat space, and this transition opens up vast new families of solutions in which the bubbles 
of these geometries change shape, or "wiggle," in ways that break the original isometries and 
depend upon internal dimensions that were previously frozen out of the dynamics. It is one of 
the primary purposes of this paper to find classes of such wiggling solutions that depend upon 
arbitrary functions. 

We therefore expect that, in the vicinity of the highly symmetric solutions counted in [HI HI], 
there will be many smooth solutions that do not preserve the isometries. Indeed, the simplest 
class of such solutions can be obtained by putting BPS supertubes of arbitrary shape [18] in 
one of the bubbling microstate solutions [19] . Through spectral flow one can then think of these 
solutions as fluctuating bubble geometries [201 EI] • These solutions depend on several continuous 
functions, and hence have an infinite-dimensional moduli space, of which the finite-dimensional 
moduli space counted in [TT| HI] is a very small, discrete subset. If one naively counts these 
fluctuating solutions, ignoring the back-reaction of the supertube^, the entropy one obtains 
from such a calculation is infinite [25]. However, if one assumes that the length of the throat 
in which the supertube sits caps off to the value expected from the dual CFT to correspond to 
the typical black hole microstates, the entropy of these supertube can potentially have the same 
growth with the charges as the entropy of a black hole [25] . 

It is also worth noting that the technique of using wigghng supertubes to create fluctuating 
bubbled geometries is a convenient mathematical device for creating an explicit, though rather 
restricted class of such fluctuating geometries. There are, of course, considerably more general 
classes of BPS fluctuating geometries, and conceivably there might even be BPS classes that 
involve fluctuating, regular surfaces with arbitrary functions of two variables or non-geometric 
backgrounds [26] . 

Returning to fluctuating geometries created by wiggling supertubes, it would be clearly in- 
teresting to use the technology developed in [TTl IH] , as well as earlier technology developed in 
[27 1 125 1 12U] to count the back- reacted wiggly supertube solutions in its inflnite-dimensional moduli 
space (as was done for two-charge supertubes in flat space in [29j). Given that the back-reaction 
of both the supertubes in flat space [I5l [161 IH] cind the ones in more complicated three-charge 
backgrounds [19] is expected to lead to smooth supergravity solution^, if the entropy one flnds 
this way is black-hole-like then this would go a long way toward establishing the validity of the 
fuzzball proposal for extremal BPS black holes. 

We have several goals in writing this paper. We begin by outlining a procedure by which one 
can construct explicit supergravity solutions that contain wiggly supertubes, and that depend 
on arbitrary functions. Even if all such solutions can be formally written down using Green 
functions [31] , constructing them explicitly is not straightforward because one needs both scalar 
and vector Green functions in multi-center Taub-NUT. The scalar ones were constructed by Page 

^Much in the same spirit as counting the D4-D0 black hole entropy in the regime of parameters where the 
DO's can be treated as probes p2l[23l [24]. 

^ These solutions are smooth only in the duality frame in which the supertubes have dipole charge corresponding 
to KK monopoles, but not in other duality frames [501 [T51 [TCI [TTl 119) . 
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for a regular multi-center Taub-NUT space [32j, but the vector ones are not known. An upcoming 
paper contains the construction of these Green functions, which are much more involved than the 
scalar ones. Even if one could write implicitly a full wiggly-supertube solution in this way, one 
cannot integrate explicitly the Green functions around an arbitrarily-shaped supertube profile. 
We therefore focus on a sub-class of wiggly supertubes where the supertube shape remains 
round, but the charge densities inside the supertube world-volume fluctuate. As we will show 
in Section [21 this class of solutions is parametrized by one arbitrary continuous function. Since 
the supertube shape is unchanged, the dipole magnetic fields of these solutions are exactly the 
same as those of the solution with a tri-holomorphic U{1) invariance. All of this makes it just 
possible to construct explicitly the fully back-reacted solution. It is thus our primary goal here 
to construct such a fully back-reacted wiggly solution. 

Constructing a wiggling supertube with one arbitrary function is obviously far less general 
than exciting all possible shape (and fermionic) modes. On the other hand, if the combined shape 
and fermionic modes of supertubes can capture the entropy of a black hole of finite horizon area, 
the number of arbitrary functions should then represent the central charge of the underlying field 
theory. Hence, one arbitrary function should be sufficient to capture the growth of the entropy 
of supergravity solutions as a function of the charges. The overall constant of proportionality 
can then be trivially restored. Hence the construction of explicit fluctuating geometries that 
depend on one arbitrary function is extremely important, because it can be used to capture the 
contribution to the entropy of the most general class of supertube fluctuations. 

A number of important mathematical and physical issues arise out of the construction of 
fully-back reacted solutions for wiggly supertubes and these lead to the other important results 
of this paper. 

First, if one wants to construct solutions that have the same charg black hole with 

a macroscopically-large horizon area and that have a hope via the entropy enhancement phe- 
nomenon to describe the entropy of this black hole, one must use ambipolar Taub-NUT or 
Gibbons-Hawking base spaces, whose signature is both -|-, -|-, -|-) and (— , — , — , — ). The five- 
dimensional solutions constructed using such apparently pathological base-spaces were shown 
in [33l El [71 to actually be smooth and horizonless. However, in order to construct wiggly 
supertube solutions one needs to know the scalar Green functions on these base spaces. 

Since ambipolar base spaces have a very significant apparent pathology, finding the physically- 
correct Green function is rather subtle. For example, one expects a Green function to diverge 
precisely when the distance between the source and the field point goes to zero but in an am- 
bipolar space the vanishing of the distance function from a given point in the (— , — , — , — ) region 
defines a co dimension-one "image" hypersurface in the +, +, +) region. However, in the fully 
back-reacted solution this point is at a finite distance from this zero-distance "image" hypersur- 
face, and hence the Green function should not diverge there. 

There are several ways in which one can try to find the proper Green function. The most 
obvious is to analytically continue the Green function on a regular Taub-NUT space first con- 
structed by Page in [35] (which we call from now on the Page Green function). As one might 
expect, such continuations are bedeviled by the "image surface" problem and, as we will show in 
Section [31 this procedure gives a Green function that solves the correct equations but diverges 
for non-coincident points. One can similarly try to use parts of the Page Green function and 
recombine them in such a way as to obtain a physical Green function, but the obvious procedures 
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do not work either. The "best" Green function that we could construct using this direct method 
has cusps at non- coincident points, and at first investigation does not appear to be good-enough 
for our purposes. It would be interesting to see whether a better-behaved Green function for 
ambipolar spaces can be constructed directly. 

The only effective strategy so far for constructing the requisite Green function on ambipolar 
spaces is to use the fact that the five- dimensional, Lorentzian geometry is regular, construct a 
Green function in five dimensions and then reduce it to the proper ambipolar Green function 
by descending to the base in four dimensions. As we will see in Section [3l this indirect method 
has the advantage of yielding the physical Green function, but has the disadvantage of being 
technically challenging for even the simplest of ambipolar spaces. We show how to use this 
method to construct the Green function for an ambipolar Gibbons-Hawking space with harmonic 
function j;^^^ — j^^^ starting from the five- dimensional Green function on global AdS^ x 5*^. We 
also compare this physical Green function to the possible extensions of the Page Green function, 
and explore the parallels and differences between the various Green functions. It is to be hoped 
that there might be a systematic procedure for adapting the Page Green function to obtain the 
physical Green function in general, but the complexity of the differences make this a challenging 
task. On the other hand, as we will describe below, we will learn from the one explicit Green 
function that we can construct that almost all of the essential physics can be extracted via careful 
use of probe solutions. 

Having obtained the Green function for the two-center, zero-charge Gibbons Hawking space, 
in Section H] we construct the fully-back-reacted supergravity solutions corresponding to wiggly 
supertubes in that space. It should also be noted that a family of wiggling supertubes carrying 
three charges was obtained in [2^ in a flat-space background. Here, by working in an ambipo- 
lar background we find three-charge solutions that form the largest known family of black-hole 
microstate solutions, parameterized by an arbitrary continuous function, and thus having an 
infinite-dimensional moduli space. We also analyze the conditions that these solutions must 
satisfy in order to be smooth and free of closed timelike curves. 

From our explicit solution for the wiggling supertube, we obtain the bubble, or integrabil- 
ity, equations that relate fluctuating charge densities and determine the location of the wiggly 
supertube. When the supertube does not wiggle, the solution can be reduced to a three-center 
solution in four dimensions, and the inter-center positions must satisfy some rather simple bubble 
equations. If the supertube changes shape, then it generically does not correspond to a point 
in the base of the Gibbons-Hawking space, and hence the regularity conditions will become 
more complicated than the simple bubble equations. Nevertheless, if the supertube is round and 
its wiggles come solely from the density modes, its location is still a point in the base, and 
one can still write bubble equations. These equations are junctionals of the density modes of 
the supertube and we will see how the ambipolar Greens functions simplify at Gibbons Hawking 
points and thus lead to relatively straightforward bubble equations on ambipolar bases. 

Having obtained these "bubble functional equations" using supergravity and our new Green 
function, we investigate, in Section the extent to which one can obtain the same physical 
data by studying the Born-Infeld action of a supertube in this space, much as one did for round 
supertubes. As a probe, the Born-Infeld dynamics of the supertube is sensitive only to the physics 
in its vicinity, and hence will only yield the bubble equation for the supertube point. However, 
one can use the fact that the bubble equations contain only pairwise interactions between various 
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centers to fish out, from tlie supertube bubble equations, tlie various parts of the bubble equations 
for the other points. This procedure actually yields all the bubble functional equations simply 
from the supertube Born-Infeld action. 

Remarkably, the two sets of functional equations agree with each other, despite being obtained 
by very different procedures: one set is obtained by going through a very technical construction 
of a physical Green function on an ambipolar space, then plugging these Green functions in the 
supergravity equations and demanding no closed timelike curves, while the other set is obtained 
by extracting and rearranging terms in the Born-Infeld action of a probe wiggly supertube. The 
fact that the bubble functional equations obtained in these two different regimes of parameters 
agree, points to the existence of a non-renormalization theorem similar to that which protects the 
bubble equations for f/(l) invariant multi-center solution^. This non-renormalization theorem 
was very useful for finding the symplectic form of and quantizing the [/(l)-invariant multi- 
center solution in jTH [T3], and we believe it will prove equally useful in quantizing the infinite- 
dimensional moduli space of the supergravity solutions we are constructing here. 

In Section E] we return to the full supergravity solution to analyze and estimate the entropy 
that can be stored in the fully back-reacted supertube. In [25] this entropy was estimated using 
the Born-Infeld action of the supertube (using the technology of [351 ES]), and was found to be 
much larger than that of supertubes in flat space. Indeed, the charges that control this entropy 
are not the electric charges of the supertubes, but some effective charges that are the sum of 
the electric charges of the supertube and of a contribution coming from the magnetic fluxes 
of the background. Hence, even a small supertube, when placed in a background with large 
magnetic fields can have a considerable entropy. The analysis in [25] involved probes and so 
did not include the back-reaction and it was found that the entropy of a supertube placed in 
scaling solution can grow arbitrarily large with increasing the depth of the scaling solution. The 
technology developed in this paper allows us to construct a fully-back-reacted example of this 
entropy enhancement mechanism, and to estimate exactly how much entropy can a supertube 
carry. We find that this entropy is to leading order independent on the electric charges of the 
supertube, and only depends on its magnetic dipole charge and on the charges and fluxes of 
the background. This shows that entropy enhancement is indeed a real physical feature of the 
complete microstate geometry. 

To summarize, in Section [2] we present our methodology for obtaining fully-back-reacted 
solutions that depend on one arbitrary continuous function. In Section [3] we explore several 
possibilities for constructing Green functions for ambipolar spaces, and find the physical Green 
function for an ambipolar Gibbons-Hawking space with two centers of opposite charge. In Section 
Hlwe use this Green function to construct the fully back-reacted solution corresponding to wiggly 
supertubes in this space, and find the functional bubble equations satisfied by the positions of the 
supertube and GH centers. In Section [5] we present a way to extract these equations from the 
Born-Infeld action of the wiggly supertube, and find that the two sets of functional equations 
agree. In Section |6] we use our results to show that the entropy of a supertube inside a bubbling 
solution depends on the charges of the background and not on its explicit charges, and thus 

■^We would like to remind the reader that these equations can be obtained both at weak effective coupling using 
quiver quantum mechanics and at strong effective coupling using the fully back-reacted supergravity solution, and 
their form is the same [13l [12] . 
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establish the existence of the entropy enhancement mechanism proposed in [25j . We conclude in 
Section 7. Appendix A contains some more details of the Green function constructed in Section 

El 



2 The general back-reacted solution with fluctuating elec- 
tric charge densities 

Three-charge solutions with four supercharges are most simply written in the M-theory duality 
frame in which the three charges are treated most symmetrically and correspond to three types 
of M2 branes wrapping three T^'s inside [311 EZ]- The metric is: 

dsl^ = - {Z^Z^Z^y-^ {dt + kf + (ZiZaZg)^ dsl 

+ (ZaZgZf 2) ^ {dxl + dxl) + [Z^Z^iZ^^) ^ {dxl} + dxl) + {Z^Z^Zf) ^ {dxl + dxl^) , (2.1) 

where ds^ is a four- dimensional hyper-Kahler metric [381 EH Ell EI] B If one imposes in addition 
an extra triholomorphic U{1) isometry, this space has to be a Gibbons-Hawking (GH) space |40j : 

dsl = V-\dtlj + Ay + Vdsl (2.2) 

with the fiat three-dimensional metric and W = V x A. 

The solution has a non-trivial three-form potential, sourced both by the M2 branes (electri- 
cally) and by the M5 dipole branes (magnetically): 

A = A^'^^ A dx5 A dxe + A^"^^ A dxr A dxg + A^^^ A dxg A dxio- (2.3) 

The magnetic contributions can be separated from the electric ones by defining the "magnetic 
field strengths:" 

e(^) ^ dA^'^ + d ( ^iL±3.] , J = 1, 2, 3. (2.4) 



Finding supersymmetric supergravity solutions for bubbles geometries reduces to solving the 
following system of BPS equations on the GH base |31j : 



□z, = iCw*4(e(^)Ae(^)), (2.5) 

dk + -k^dk = ZjQ^ . 

In these equations, *4 is the Hodge dual in the four- dimensional base space and Cjjk = I^ijkI- 
The operator, □, is the Laplacian in the four- dimensional metric. 



''This metric can have regions of signature +4 and signature —4 [331 IHl [3 [31], and for this reason we usually 
refer to it as ambipolar. 
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Our purpose is to solve this in a similar manner to that described in [31j , except that we want 
to allow varying electric charge densities. We will also later focus on three-charge geometries 
containing supertubes and so we will work in the D1-D5-P duality frame (where supertubes 
source smooth solutions). In this frame, the metric takes the form: 

dsj,^ = -—^=={dt + kf + ^z;z;dsl + ^={dz + A^'^f (2.6) 

^3V^1^2 V^1^2 




+ \ — {dxi + dxQ + dxj + dxl) , (2.7) 

V Z2 

where ds\ is still the metric of the GH base. In the M-theory frame, the electromagnetic fields 
all appear in the 3-form Maxwell potential f l2.3p while in the IIB D1-D5-P frame one of these 
fields has become part of the metric (12.71) . 

The first step is to solve the first BPS equation in exactly the same manner as for GH base 
spaces, [m m [71 [2] and for this it is convenient to introduce the vielbeins: 

= V-^(# + A), e'^+i = v'^dy\ a = 1,2, 3, (2.8) 

and the two-forms: 

= A 6"+^ ± i tabc e^^^ A e'+^ , a = 1, 2, 3 . (2.9) 

One then takes 

e(^) = -f2{da{v~'K^)) (2.10) 

a=l 

for some harmonic functions, , on the M.^ base of the GH space. We specifically require the 
magnetic fluxes to be independent of ip. These field strengths have potentials: 

B'^^'^ = V-^K^ [dij + A) + i^^^ -dy, Vx^"^^)= -VK\ (2.11) 

where O^^) = dB^^\ 

The source of the second BPS equation is independent of ip and so the inhomogeneous solu- 
tions for the functions Zj follow the standard form: 

Zi = \CukV-^K'K'' + Li, (2.12) 

and the functions Lj are still required to be harmonic but are now going to be allowed to depend 
upon all variables including the fiber. Thus we have 

□Lj = 0. (2.13) 

As before, there is a natural Ansatz for the angular momentum one-form k: 

k = ^{di{j + A) + u ■ dy . (2.14) 



8 



The third BPS equation then reduces to: 

{fiVV - VVfi) + V[d^uj + Ad^fx] + V x uj + {d^uj x A) 



^ Z/V(V-^i^^) , (2.15) 



1=1 

where A is the vector field in the fibration of the GH metric and satisfies V x A = W. 
If one defines the covariant derivative: 

V = V - Ad^, (2.16) 

then one can write fl2.15p as 

3 

(/iW - VVjj,) + Vxuj + Vd^oo = -VJ2 ZiV{V-^K^) . (2.17) 

1=1 

Moreover, the four- dimensional Laplacian can be written: 

□F = V-^[V'^dlF + V-VF]. (2.18) 

The third BPS equation has a gauge invariance: k ^ k + df and this reduces to: 

/i^/i + (9^/, u ^ u + Vf , (2.19) 

The Lorentz gauge-fixing condition, d *4 /c = 0, reduces to 

V^d^fx + V-uj = 0, (2.20) 

and we will adopt this gauge throughout. 

Now take the covariant divergence, using V, of (I2.17p and use the Lorentz gauge choice, and 
one obtains: 

3 



V' an = V ■ (v J2 ^I'^i^'^^^)) ■ (2-21) 
1=1 



Remarkably enough, this equation is still solved by: 

= IV-^CijkK^K-^K^ + \ V-^K^Li + M, (2.22) 

where, once again, M is a harmonic function in four dimensions. Finally, we can use this solution 
back in fl2.17p to simplify the right-hand side to obtain: 

Vxu + Vd^uj = VVM - MVV + ^{K^VLi - LiVK^). (2.23) 

Once again one sees the emergence of the familiar symplectic form on the right-hand side of this 
equation. One can also verify that the covariant divergence (using P) generates an identity that 
is trivially satisfied as a consequence of VV = V x A, fl2.20p . f l2.22p and 

nil = DM = 0. (2.24) 

While it might be possible to find, we do not have an explicit closed form for uj, but we will not 
need it. 

We have thus generalized the results of [39l HH [6l [7] to completely general, fiuctuating electric 
charge densities. 
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3 Scalar Green functions on a GH space 



In order to solve f l2.24p for varying charge densities, one needs the scalar Green function on 
a GH space. When the GH space is a regular Euclidean space (the function V is everywhere 
positive), the Green function is well-known [32]. However, as we explained in the Introduction, 
when the GH space is ambipolar finding the Green function is not so straightforward, and the 
simple analytic continuation of the Green function on a regular GH space does not produce a 
physical ambipolar Green function. We therefore turn to an indirect method for constructing 
this Green function: we first construct the Green function for a five-dimensional solution that 
has an ambipolar base space: AdS^ x S'^, and then we reduce this Green function. 

This section reviews some of the known material about Green functions on GH spaces and 
the relationship between AdS^ x 5*^ and an ambipolar two-centered GH space. Much of the 
rest of this section contains the (highly technical) derivation of the Green function that we will 
subsequently use. Those who wish to skip most the technical details can get a good idea of the 
general strategy from Section 13.31 The key physical results that we will need about the Green 
function are contained in Section 13.71 which discusses the behavior of the Green function in 
various limits. In particular. Section 13.7.21 contains the essential result that demonstrates that 
the bubble equations at the GH fibers (but not at the supertube itself) are unchanged by the 
fluctuations and only depend upon the total charge of the fluctuating component. 



3.1 Review of known Green functions on regular GH spaces 

The Green function on a positive definite, multi-center Taub-NUT space have been known for 
some time. Consider the metric (12.21) with V harmonic and everywhere positive: 



V 



qn 

I 



with Qn > 0. The Green function was found by Page in [32j, and is given by: 

1 sinh U 



where 
and 



IGvr^A cosh U — cosT 
A = \x — x'\ 



T 
U 



arctan 



cos 



tan 



n 



COS 



ij^g^logj 



(3.1) 



(3.2) 
(3.3) 



(3.4) 



and for brevity we have introduced spherical polar coordinates (r^. On-, (pn) around each pole at 
position Xn- The only singularity of this Green function is at ?7 = T = 0, which for regular GH 
spaces with Qn > only happens when the points are coincident: ip = ip' and x = x'. 
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By integrating this Green function against the fiber ijj, one obtains the trivial Green function 



on M^: 



47r 







dtljG = ——. (3.5) 
^ 47rA ^ ^ 



One might hope that that equation (13.21) would also give the Green function for ambipolar 
spaces where some of the g„ are negative, and indeed the Green function one obtains (which we 
refer to as the Page Green function) satisfies the appropriate differential equation. As before, 
this Green function is singular when U = T = 0, but in an ambipolar space this does not only 
happen when the two points are coincident. Indeed, as one can see from Eq fl3.4p . for a fixed 
{ip'jx'), the points ('?/', x) defined hj U = T = belong to a co dimension-two hypersurface, and 
in general U and T can be both positive and negative. 

One could also calculate the integral of the ambipolar Page Green function over the fiber ip: 

dip G Page = TTTT-r—r (3.6) 



|?7|47rA 

and see that this integral is discontinuous along the surface U = 0. 

One can also try to construct a physical Green function on an ambipolar space by piecing 
the Page Green functions on the patches U > and U < 0: 

G th= ^ sinh|?7| 

patch Xg7^2^ cosh U — COS T 

The ifj integral of Gpatch then gives the correct three-dimensional propagator, but this function 
has a non-physical cusp across the [7 = hypersurface. Apart from Gpage and Gpatch, there 
does not seem to be any natural guess for a physical continuation of the multi-center Taub-NUT 
Green function to ambi-polar spaces, and hence we do not have a direct method to construct 
this Green function. 

Having failed in the brute-force construction, we now take the less direct road by remembering 
the ultimate goal of the overall construction: A smooth five-dimensional geometry. We there- 
fore construct the Green function for a particular family of ambipolar spaces by dimensionally 
reducing the Green function of the smooth five-dimensional solutions built using these spaces as 
a base. 



3.2 Preliminaries: the simplest metric with an ambi-polar base, 
global AdSs x 

It has been noted by several authors (see, for example, |23lll2]) that the five-dimensional metric 
metric arising from an ambi-polar Eguchi-Hanson metric with GH charges +q and —q is, in fact, 
a 1>q quotient of global AdS^ x S^. For simplicity, we will locate these GH charges on the z-axis 
at z = ±a and define: 

r± = + izTa)\ (3.8) 
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where {z, p, 0) are cylindrical polar coordinates on the M.^ base. The harmonic functions for this 
solution are 

V = q(- - -), K = k(- + -), (3.9) 

\r-i_ r_/ Vr_|_ r_/ 

L = , M = + --(- + - , (3.10) 

where the constant in M has been chosen so as to make the metric regular at infinity. 
The vector potentials for this solution are then: 

V r+ r_ J aq r^r^ 

The five- dimensional metric is then: 

dsl = - Z-\dt + fi{dij + A) + + Z{V-\d^ + A)^ + V{dp^ + p^d(f)^ + dz^)) , (3.12) 
where 

4 k^ 1 

Z = V-'K^ + L = (3.13) 

q [r^ — r_) 

p = V-V + IV-'KL + M = ^ + - ^. 

g"^ (r+ — r_)"^ a 

To map this onto a more standard form of AdS^ x 5*^ one must make a transformation to 
oblate spheroidal coordinates like the one employed in [13] to map positive-definite two-centered 
GH space onto the Eguchi-Hanson form: 

z = a cosh2^ cos^, p = a sinh2^ sin^ , ^>0, < 6 < tt . (3.14) 

In particular, one has r± = a(cosh2,^ =p cos 6*). (Note that we have taken the argument of the 
hyperbolic functions to be 2,^ for later convenience.) One then rescales and shifts the remaining 
variables according to: 

r ^ i§^t, V^i^^V^-l^t, + (3.15) 

and the five- dimensional metric takes the standard AdS^ x S"^ form: 

dsl = Rl[- cosh^ ^ dr^ + di'^ + sinh^ ^ d^f^ + Rl [dO'^ + sin^ 6 d^l] , (3. 16) 

with 

= 2R2 = Ak. (3.17) 

Note that the metric involves global AdS^ with — 00 < r < 00. 

One should also recall that the GH fiber coordinate has period An and therefore, for |g| = 1, 
the angles, ipj, both have periods 27i. For g 7^ 1, the Z\q\ orbifold associated with the GH points 
emerges as a simultaneous Z|q| quotient on the longitudes of the AdS^ and S"^. 



12 



It is instructive to observe that a particle that is at a fixed spatial point in the GH base has 
a world-line with constant ^ and 6 and 



dr ' dr 



dsi, = - 16k' cos^ e dr'. (3.18) 



This is the world-line of a observer rotating around the AdS^ and 5*^ simultaneously along a 
curve that is null for ^ = | and otherwise time-like. The surface ^ = | is the simply the critical 
surface defined by = in the GH base. As has been noted elsewhere, and is evident from 
( I3.16p . the five- dimensional metric is completely regular across this surface and the only artifact 
of an apparently singular four-dimensional base is that a stationary particle on this base is simply 
following a null curve on the critical surface in the five- dimensional geometry. 

3.3 The strategy: reducing a Green function from five to four dimen- 
sions 

Our goal here is to obtain the scalar Green function for the ambi-polar GH base considered above, 
and construct this Green function in such a manner that it can be used to create five- dimensional 
solutions by the methods outlined in [31j . One can obtain such a Green function from the Green 
function for AdS^ x 5"^ by integrating it along the time coordinate, t, in f l3.12p . 
To be specific, suppose that we have a function, Q{x]x'), that satisfies 

a,g{x;x') = -^S{x;x'), (3.19) 

on AdSs X S'^. Let x = {y,t), where y stands for the coordinates on the GH base. Write the 
metric fl3.12p as: 

dsl = - Z-^{dt + k^dy"'Y + Zhmndy'^dy'', (3.20) 
and observe that f l3.19p is equivalent to: 

-Z-'d^g{x-x') + l={d^ - k^dt) {Vhh^^idn - kndt)g{x;x')) = -^6{x-x'). (3.21) 

Since all the metric coefficients are time-independent, one can take all the time derivatives 
to be total derivatives of the various terms in fl3.2ip . Integrating this identity over t then yields 
the identity of the form: 

ayGiy;y') = -^d^ {Vhh^^d^Giy;y')) = ^ , (3-22) 

provided that one can drop the boundary terms from the t integration. We will construct the 
Green function with a periodic time variable and then analytically continue using t — t' = iX. 
With this prescription the boundary terms will indeed vanish. We therefore define 

/oo 
dXgiy,T = T' + tX;y',T'), (3.23) 
■oo 
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and this will give us the required Green function on the GH base. One should note that because 
of the off-diagonal "angular momentum" terms in the metric (13.201) . the analytic continuation of 
G{x;x') will not produce a real integrand in (13.23^ . however, the imaginary parts will disappear 
as boundary terms in the integral. Periodic identification of the time coordinate, r, introduces 
an infinite number of image sources into the original AdS^ space and one might be concerned 
that this will generate spurious sources in G{y]y'). However this is not a problem because we 
are integrating over the periodic variable and the periodicity is commensurate with that of the 



3.4 The scalar Green function on AdS^ x S'^ 

The scalar Green functions for AdSs x S"^ has been extensively discussed in [S] and here we 
simply follow their prescription. It also turns out that the particular AdS^ x S"^ given by (I3.17p 
has some extremely nice properties with respect to the procedures of [Hj. 
Consider the more general equation: 

{DAds, + Os^-M^)g{x;x') = -^6{x;x'), (3.24) 

where the factor of i has been included so that, upon passing to the Euclidean version, the Green 
function is real. Expand the Green functions into spherical harmonics on S"^ one obtains: 

g{x; x) = ^ V QiA^; w') Y,,UO, V2)Yl^{e\ ^'^) , (3.25) 

where w and w' are coordinates on AdS^ and 

(oA,,^-M'-^-^^±^]g,^Uw;w') = — ^ 6{w-w'). (3.26) 

Since the functions, Qi^m, do not, in fact, depend upon m, we may perform the sum over m 
to obtain: 

G{x;x') = ^J2g,{w;w'){2e+l)Pi{cos^), (3.27) 

where 7 is the great-circle angular separation on the S"^: 

cos 7 = cos 9 cos 9' + sin 9 sin 9' cos((y92 — (p2) ■ (3.28) 

One can now solve the equation fl3.26p for in terms of the "chordal distance" on AdS^. 
Consider the hyperboloid: 

[Xof + {XsY - {X,f - {X^f = Rl , (3.29) 
with the metric induced from 

ds^ = - dXl - dXl + dXl + dXl . (3.30) 
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With the coordinate change: 

Xq = i?icosh^ COST, X3 = i?icosh^ sinr, (3.31) 
Xi = i?isinh^ cosy^i, X2 = -Risinh^ siny^i, (3.32) 

the induced metric becomes precisely the AdS^ factor of fl3.16p except that here the time coordi- 
nate, r, is periodically identified with period 2n. As we noted above, this periodic identification 
is consistent with the the world-lines fl3.18p over which we wish to integrate. 
The chordal distance is given by: 

y = _(Xo-X;)2-(X3-X^)2 + (Xi-X()2 + (X2-X^)2 = 2Rl{C - 1), (3.33) 

where 

( = cosh ^ cosh ^' cos(r — r') — sinh^sinh^'cos((y9i — y^'^) . (3.34) 
One can then show that [H]: 

where F is a hypergeometric function and 



A = l + ^J^+i{i + l) + M^Rj. (3.36) 
K2 y Ki 

We want the Green function for M = and Ri = 2R2 for which there is a remarkable 
simplification: A = 2{i + 1) and the sum in ( 13.27^ can be performed explicitly. Let rj = 
(C + a/C^ ~ 1) = C ~ a/C^ — 1 and observe that f l3.27p has the form: 



1 rj 



Sn^RiRl 1 - ri^ 

7f d 



^(2£ + l)r72^+ip,(cos7) 



/ ^ 2/4-1 / \ 



£=0 

n2 



— (77(1 + 77 —2ri cos7j 



M'^RiRl l-T]^ dr] 

1 iv + v'^) 



We therefore find 



M^RiRl (r72 + r/-2-2cos7)i 
1 C 



8y2 7r2i?ii?2 (2C2_ (1 + COS7)) 



3 1 
2 



(3.37) 



(3.38) 



where 7 and C, are given by fl3.28p and (13.340 . The result is thus a relatively simple combination 
of chordal distances on AdSci and S"^. 
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One should note that in the Euchdeanized AdSs with r — > ir, it follows from fl3.34p that 
C > 1 with equality if and only if ^ = r = r' and = fi- The denominator of f l3.38p only 
vanishes if = 1 and 7 = 0, which means that only singularity of Q occurs when the points 
coincide on the Euchdeanized background. In addition, C > 1 means that rj < 1 and this is the 
condition for convergence of the series in (13.371) . 

We have therefore constructed the scalar Green function on Lorentzian AdS^ x S"^ with peri- 
odic time and seen that it has the appropriate singular structure when Wick rotated to Euclidean 
AdS^ X S"^. To get to the Green function on the GH base we have to use the transformation 
(I3.15P and this will lead to a somewhat different analytic continuation. We now investigate this 
in more detail. 

3.5 Reducing the scalar Green function to the GH base 

To get the requisite Green function, G{y;y'), on the GH base, we have to make the changes of 
variable: 

if, = ^^^~T, if2 = 0-^^ + 2r, (3.39) 

make the analytic continuation t — t' = iX and perform the integral (13.231) . Simplifying the 
normalization for the present, we have: 

/'OO /"OO /■ 

X = 8V27r^RiRl / G dX = / ^ j dX , (3.40) 

(2^^ — (1 + cos 7)) 2 



00 



ai cosh A + i6i sinh A /„ 

■dX, (3.41^ 



(02 cosh 2A + ib2 sinh 2A + oq) 



where 



where the and bj are real variables. It turns out that this integral is elementary and it is given 
by: 

X = ^ + / ^ (3.42) 

sj\{a2 + ib2) y|(a2-«&2) 

p = I (oq - (02 + 62))"^ [(«o - a2)«i - + ^((«o + «2)&i - ai&2)] , (3.43) 

and p is its complex conjugate. The result, while apparently rather simple, is an extremely 
compressed form of the answer we are seeking and we first rewrite it in terms of the coordinates 
of the points. 

Define the following combinations of coordinates: 

u = cosh ^ cosh ^' , V = sinh ^ sinh ^' , x = cos 6 cos 6' , 

y = sin sin 0', z = e^^''^~^'\ x = (0-0'), (3-44) 

then one has 

fli = u — ^v{z + z~^), bi = ^v{z — z~^), 

ao = al + bl — {l + x), 02 + 262 = {u — vz^ — ye~'^^z . (3.45) 
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It is also convenient to take: 



A = z{al- {al + hl)) , m = 2{l + x)uv + yiyU^e''^ + v'^e'^) 

V = (cto - ct2)ai - + + 0-2)^ - aife2) , 

V = (2(u-v?- 



(3.46) 
(3.47) 

1 + x + y)) {2{u + vf -{1 + x- y)) + 8uvy{l - cos x) , (3.48) 

V] . (3.49) 



z± = ^[y'^ - (1 +xf + 2{u^ + v^){l + x) + Auvy cos x± a/(1 + x)2 - y 
Adjusting the overall normalization, the desired Green function is then: 



G = 



I TT 



2-R1-R2 



gdx 



Re 



z V 



a Aa/oT+TS^J 



(3.50) 



This is the Green function on the GH base and our task is now to unpack the details and 
understand the result. As with all Green functions, the physics lies in the singular structure and 
asymptotic behaviors. 



3.6 Properties of the Green function 

While the Green function fl3.50p is rather unintuitively expressed and represents a very compli- 
cated function of the physical variables, its overall structure can be understood fairly easily. Here 
we will summarize the important features and relegate some more technical details to Appendix 
A. Since z is simply a phase, the physical domain has \z\ = 1, however because we are going to 
want to use residue calculus to evaluate integrals over physical values of 2, we are going to want 
to understand the behavior of this Green function in the domain \z\ < 1. 

Since all the constituent functions are analytic, the only potential sources of singularities are A 
and \/a2 + ib2- The bottom line will be that the function A has essentially the same singularities 
as the Page Green function, while the remaining components play the crucial role of canceling 
the unphysical singularities or cusps at the critical (V = 0) surface and at the dangerous "image" 
surfaces, and giving a smooth physical Green function. 



3.6.1 Returning to the GH geometry 

We start by recasting our result, as far as we can, in terms of the original GH geometry. Following 
[32], we define r±, r±, 9± and 9'^ to be the radial distances and polar angles of the two points as 
measured from the GH points aX z = ±a, and we define the functions: 

T = + ^iarctan 

± 

where 

- ^\ = (r| + - 2 r±r^ (cos ^± cos + sin ^± sin ^^cos(0-0')))^- (3-53) 



cos 2! 



COS 



-e' 



tan ■ 



(3.51) 



r-t + ri|_ + A 
r-t + r j_ — A 



(3.52) 



A = f 
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Note that these functions are normahzed differently from those of [32j in that we have divided 
by the GH charge, q. 

It is a rather tedious process to relate the foregoing, more geometric, functions and variables 
to those of the previous subsection. The basic relationship is that r± = a (cosh 2^ =p cos^) and 

= a(cosh2^' =p cos^'). Then by drawing triangles and finding expressions for various side 
lengths one can obtain: 

^r+r_ cos I (6'+ — 6_) = a sinh2^ , ^yr+r. sin |(^+ — 6'_) = a sin ^ , (3.54) 

^Tj^t^ cos \ {Q^ + 6'_) = a sinh 2^ cos Q , ^r+r_ sin 1(6*+ + 6'_) = a cosh 2^ sin 6' . (3.55) 

The function. A, is manifestly a quadratic in and one can easily show that it has the 
following form 

A = miz — Zj^Mz — zJ) with Zj^z^ = — . (3.56) 

m 

One can then establish the following identities: 

e^^ = ^z, = ^z^, A = aVV, 

\m\ \m\ 

1_ [(r+ + r'^f - A^] ^ [(r_ + T'_f - A^] ^ (3.57) 



From which one obtains: 



4a 



sinht/ = — A (cos ^ + cos ^') , (3-58) 

2 a \m\ 



and 

1 1 a 1 sinh U 



A 2 |m| (cosh f/ — cosT) A (cos 6' + cos 6'') (cosh [/ — cos T) 



(3.59) 



Thus we have recovered the core part of the Page Green function. 

We have not found any particularly simple form for 02 + 0)2 or v in terms of U and T, but 
as we will now see, these functions play an essential role in avoiding jump discontinuities across 
critical and image surfaces. 

3.6.2 Residues and pole structure 

It follows from f l3.56p that \zj^z^ \ = 1 and so if z± is inside the unit circle, then z=p is outside the 
circle. One can also verify that 

z/ = — {1 + x){u — vz) — y{v — uz)e~^^ . (3.60) 

and 

(02 + ^62) = {{l + xf -y^y\u^ + ye-'''A) = {cosO + cosffy^^u^ + y e-'^'k) . (3.61) 

The importance of this identity is to show that there is a uniform, analytic choice of square root 
on the locus A = 0. More to the point, the roots A = correspond to: 

z = z± V = ± (cos 9 + cos 9') \J a2 + i&2 • (3.62) 
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Thus, we have: 

(cos + cos ^^') ^3^3^ 



V 



z=z± a m [z^ — z. 



Res 

\-aAy/a2 + ib2 

Note the absence of ± signs in the result. 

The factor of ^^^'^^^ thus has two very important effects on the result. First, it cancels the 
(cos^ + cos 6'') in (I3.59P to give the Page Green function. If this term were not canceled then 
there would be the spurious "image" poles, alluded to earlier, at 9 = n — 9'. The second effect 
is more subtle: there is an extremely important effect on signs coming from (I3.62p . Indeed, the 
effect is precisely the difference between: 

± — , — /t T- = — , (3-64) 



27ri Jc {z - Zo)iz + zo) 2zq 2m Jc {z - zo){z + zq) zq 2zq 

where the sign depends whether the (counterclockwise) contour, C, surrounds ±2:0. The sign 
of the second integral does not depend upon whether the contour surrounds either +zq or —zq. 
This removes jump discontinuities in Green function. 

3.7 Solutions and their asymptotic behavior 

In this section we set up our conventions for constructing the wiggling solutions from source 
currents and then examine the asymptotic behavior of these solutions in various limits. These 
limits provide the crucial physical data for the regularity conditions on the wiggling supertube 
and in the bubble equations. 

3.7.1 Short-distance behaviour 

First, in the limit that f f' and ip ^ ip' one has, from (13.521) . T, [/ — t- and 
z a cos 9 sinh U p 2 



A A (cosh [/- cos T) ' y/aPi^ y/^^ 

Thus 



(3.65) 



- 1 cos^ sinhf/ 1 |sinht/| 2 \U\ 

G ~ ~ ~ . (3.66) 

A I cos 0| (cosh U - cos T) A (cosh U - cos T) A {U^ + T^) ^ ' 

The absolute values come from the a/ cos^ 9 and have the effect of keeping the coefficient of the 
singular part of the Green function positive as one crosses the surface ^ = 0. 

For a small separation in the GH base, the infinitesimal proper length is given by: 

(Ss)^ = 1/-^((V^'-^) + g(cos^+-cos^_)(0-0'))^ + VA. (3.67) 
One then finds that f l3.66p can be written as 

S ~ (3.68) 
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Note the absolute values in the denominator. As we will discuss below, this result means that G 
is the properly normalized Green function for our purposes. 

The Green function, G, depends upon ip, ip' via cosT, where T is given by (13.52^ and so G is 
periodic on the interval [0,47rg]. It is, of course, elementary to produce a Green function with a 
reduced periodicity, like Arc, simply by summing over image sources appropriately. Here we will 
work with the function G and sources that have the same periodicity and so the solution to the 
Laplacian with a source density, p{ip'), along the GH fiber located at r ' is therefore given by: 

$(f,^;f')= / G(f,^; f',^')p(V'') (3-69) 







It is convenient to define the associated total charge by: 

Q = pi^') dijj'. (3.70) 

Jo 

In the coincidence limit, when one approaches the source, r — )■ r' ^ ip' the dominant contri- 
bution to f l3.69p comes from integrating through the singularity and so 

$ - / ^- .php') #' 

V, |r-i(^-^')' + V^\^^^' 

This is precisely what one should expect: For a constant charge, the charge density is p = ^ 

and thus (13.711) gives $ ~ ^, which is the canonical form in the GH base. Moreover, if one 
approaches any line distribution of charge then is should appear, locally, like a constant charge 
and hence one must obtain (13.711) . This is the reason for selecting the normalization of G in 



3.7.2 Asymptotics at the Gibbons-Hawking points 

Another important limit is to understand the behavior of G and the solution $, defined in 
(I3.69p . near the GH points. This evaluation is straightforward and one finds that at r± = 0, or 
r = (0, 0, ±a), one has 

G{f={0,0,±a),^-f',^') = ^. (3.72) 

Observe that this is independent of ip and and so the solution (I3.69P depends only on the 
total charge of the supertube, and not on the way this charge is distributed on the supertube 
world- volume: 

$(f=(0,0,±a),V';r',V'') = - / P(V^')#' = ^- (3.73) 

r± Jo 

One should have anticipated a result of the form (I3.73P from the outset: At the GH points 
the GH fiber, parametrized by ip collapses to a point but the space-time looks locally smooth. 
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Regularity therefore requires that any non-trivial Fourier mode on the fiber must vanish at the 
GH point and only the constant mode can survive. Thus the GH points only feel the total charge 
of the fiuctuating solution. 

Now recall that the bubble equations at the GH points are local conditions and so it follows 
that in a geometry with fiuctuating electric charge densities along the GH fiber, the bubble 
equations for the GH points are unchanged by the fiuctuations and only depend upon the total 
charge of the fiuctuating component. 



3.7.3 Behaviour at infinity 

We now consider the behavior of the Green functions and solutions of the form (I3.69P as r = |r| 
becomes large. This corresponds to taking ^ to be large and (13.141) implies 

e-2« ~ — . (3.74) 

One can now take the limit of large ^ in (I3.50p but it is simpler to take the large ^ limit directly 
in (I3.4ip . In particular, note that in the large ^ limit, with the analytic continuation r — r' = iX, 
one has: 

C ~ I e^(cosh,^' cos(r — r') — sinh^' cos((y9i — 99'J) (3.75) 
ie^(aicoshA + ibisinhX). (3.76) 
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Thus 



G = / ^ jdX (3.77) 

y_oo (2(2- (l + cos7))i 

00 -j^ poo 2 

d^ = Tz 7:—::r^--r—rT^dX (3.78) 



4 a 4 a J _^ {ai cosh X + ibi sinh X)"^ 

1 2e-2« 



2 a {ai — ibi){ai + ibi) a (cosh^' — zsinh^')(cosh^' — z ^sinh^') 
1 z 



r sinh ^' cosh ^'{z — coth ^') {z — tanh ^') 
Now suppose that a charge density has the Fourier series: 



(3.79) 
(3.80) 



p{^P') = (3.81) 

n=— 00 

for which reality implies a_„ = a^. Then one finds that the asymptotic form of the fundamental 
solution, f l3.69p . is given by: 
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This integral is an elementary residue for n non-positive, while the result for n positive is the 
complex conjugate of the result for n negative. The solution is thus: 



+ ^2 i'^nW"' + anW'^) 



n=l 



Observe that we have the following limits 



$(o,v) 



w 



tanhC'e*'^. (3.83) 



(3.84) 



and ^' = puts the supertube on the axis between the two GH points, while ^' — )■ oo corresponds 
to moving the supertube off to infinity in the GH space. Thus the asymptotics of the solution 
averages out if the supertube is between the two GH points and is directly proportional to the 
density if the supertube is at infinity. 

We are generically interested in putting bounds upon asymptotic values of our solutions, 
and to that end we note that the asymptotic functional dependence fl3.83p satisfies the Laplace 
equation on the semi-infinite cylinder defined by — log(tanh(^')) and ip. It follows from the 
maximum modulus theorem that the asymptotic behaviour of the fundamental solution as — ?■ oo 
is bounded by the boundary values, f l3.84l) . 



4 The fully back-reacted solution for a wiggling super- 
tube 

It now remains to assemble all the constituents of the fully back-reacted solution and to check 
the regularity at the supertube and at the GH points and then impose any constraints required 
by regularity and absence of CTC's. We found, in Section [3. 7.2[ that in a geometry with fluctu- 
ating electric charge densities along the GH fibers, the bubble equations are unchanged by the 
fluctuations and only depend upon the total charge of the fluctuating component. This has the 
important consequence that regularity at the GH points for the corresponding wiggling supertube 
has exactly the same form as the regularity conditions for the round, non-wiggling supertube. 
There are, however, two important differences. 

The first difference is that around the wiggling supertube the standard regularity conditions 
are replaced by a simple generalization to a local regularity conditions along the fiber. Sec- 
ondly, we find that the dipole-dipole interactions will allow large negative values of the angular 
momentum, J, of the supertube, and thus enable entropy enhancement, but this also has the 
concomitant effect of having the supertube approach the critical {V = 0) surface where it be- 
comes larger and floppier, while remaining round. The fluctuating charge density thus generates 
a displacement in the base space and an associated change of the supertube radius. We can then 
interpret the result as a shape mode that preserves roundness: The fluctuations in charge density 
can be thought of as being the result of differential expansion of different parts of the supertube 
(while keeping it round) diluting charge in some regions and concentrating it in others. This 
process is depicted in Fig. HI 

We now construct the solution in detail and exhibit all the regularity constraints. 
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Figure 1: A round supertube can retain its round shape but concentrate or dilute charges in different regions as 
it expands or contracts. Some large sectors of the original tube can be mapped into smaller sectors, concentrating 
charge, while some sectors can become expanded resulting in a charge dilution. Even though the shape of the 
supertube remains circular, the charge density modes induce a radius change and can thus be seen as special class 
of shape mode. 



4.1 The constituents of the solution 

As shown in Section [21 the full solution for a round supertube with fluctuating charge densities 
inside a GH space is given by 8 harmonic functions in five dimension: 



r+ r_/ Vr_|_ r_/ A 

2 



+ ^ A=l 

M = mo + m(— + — ) + j{f,tp;r'), 



where the components of the solution sourced by the fluctuating densities are 

A^(f, ^; r') = / G(f, r', ^') d^^' , 

Jo 

j(f, V-; r ') = / ' G(f, ij- f xlj') p{iP') d^P' , (4.2) 
Jo 

and ^ ^ 

= - TT^iJKkjkK , m = ——^CijKkikjkK ■ (4.3) 

Note that in (14.11) we have chosen a gauge so that the have equal charges at the GH points. 
For simplicity, we will also assume that g > so that integrals like (14. 2 p run in the canonical 
(positive) orientation. We will also assume that the kj are all positive so that it is simple to 
assign signs to square roots in the metric and other functions. 
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With these choices we have: 

Z. = + + ,4.4) 

^2 = 7 r + -r + Aa, ^3 = 7 r, 4.5 

= mo + j + 5 + ^— -(F - i) + — (fciAi + /caAajF, (4.6) 



where 



F^.il±±lA (4.7) 

r+ — r_ 



It is also convenient to define a special value, F, of F at the supertube: 

F ^ F,^, = - P^ll^ , (4.8) 
(r^ — r_) 

where r'^ = — f±\. 

Note that we do not have an explicit solution to the (linear) equations, fl2.20p and f l2.23p . for 
oj, but we will not need the explicit solution in the following. 



4.2 Supertube regularity 

Consider a supertube located at r = r ' in the base of the GH space and with charges 
associated with K^, Li, L2 and M. That is, suppose that only these functions diverge as A~^ 
for A = \ f— f'\ — )■ 0. Moreover, as above, we allow Li, L2 and M to have densities that depend 
upon -0. The regularity of this supertube amounts to requiring that the metric (12. 7p be regular 
and without CTC's at the supertube. There are two potential dangers: (i) Divergences parallel 
to the fiber, and (ii) Dirac strings in u. As discussed in for this supertube, the former 
condition requires that for all ip: 

lim A'^lZsiK^y - 2ijVK^ + Z1Z2V] = 0, (4.9) 

A— !>0 

where A = \ f — f'\. 

To analyze Dirac strings, consider any two-sphere, S*^, in the base of the GH space. If 
the magnetic field, V x cJ, has a monopole source inside the S"^ then u will have a Dirac string 
somewhere on this S'^. Let (p be the polar angle around this Dirac string and suppose that the 
charge density is independent of then one can let ip vary in such a way that (dip+A) vanishes as 
goes around the string. If one does this, the closed curve around the Dirac string is necessarily 
time-like in the metric (12. 7p . If the charge density varies then the foregoing argument may not 
work, but if the variation is small then there will still be CTC's associated with the Dirac string. 
It is conceivable that one might find a conspiracy for large variation in the charge density and 
that one might avoid CTC's associated with Dirac strings, but it is simpler to require that there 
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are no Dirac strings in lo no matter how it depends upon That is, for all '\\) we impose the 
condition: 

/ (V X - n = 0, (4.10) 

where the S'^ is any small sphere in the base drawn around the charge source and n is the 
unit outward normal to this S"^. 

Suppose that the charge density at r = r ' has a monopolar magnetic field and consider a 
vanishingly small S*^ around r' . As A — )■ 0, the Cartesian M'^ components of Q will diverge, in 
a direction-dependent manner, as A~^. Also assume that r' is not a GH point and so V and 
A = A ■ dy are smooth at r = r '. Now imagine integrating equation fl2.17p over this 5*^: The 
terms involving d^Q and A x d^Q and VAd^fi do not contribute to this integral as A — > 0. 
Therefore we find that 

f (Vxa;)-n = f V (Vij. - ^ V(V-^K^)) ■ n . (4.11) 

For this to vanish as A — )■ one must have: 

hinA[Vn - Z^K^] = 0. (4.12) 

This result is of the same form as the result for constant charge densities, but here this equation 
is to be applied as a functional constraint that must be true for all ip, just like f l4.9p . 
One can use f l4.12p in fl4.9p to obtain the simpler condition: 

lim A^TV'ZiZs - ZsiK^?] = 0. (4.13) 

It is slightly more convenient to work with f l4.12p in f l4.13p . 

Let pi{ip), Pii'ip) and p{ip) denote the charge densities that generate the functions Li, L2 and 
M via the expression fl3.69p . Let k be the magnetic dipole charge of the supertube, that is: 

~ A, (4.14) 

as A — > 0. In addition, let V, K^, K^, L3 and denote the (finite, ■^/'-independent) values of 
the corresponding functions at the supertube location {-r = f'). Using the limiting behavior in 
(13. 71 p . one finds that (14.121) becomes 

47rg[irVi(V^) + K^P2W + 2 V'p(V^) ] = kL^. (4.15) 

Similarly, f l4.13p becomes 

[Airgpiitfj) + kV"^K^] [ATTqp2{ip) + kV-^K^] = Pv^^Za . (4.16) 

Combining this with ( I4.15P gives the simple relationship between the densities: 

kp{^) = 27rgpi(^)p2(V^). (4.17) 
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Note that if one specifies pi{ip) then these identities fix P2(V^) and p{ip) in terms of ip- 
independent geometric factors. 
Finally, define the charges: 



Qa = 



J 



(4.18) 



It is useful to note that (14.151) integrates to: 

[K^Qi + k^Q2 + 2VJ] = kLs. (4.19) 
For constant charge densities one thus gets this equation combined with (14.171) . which yields 

Qi Qi 



J 



2k 



(4.20) 



4.3 The bubble equations 

The bubble equations at the two GH points are given by setting /i = at each point: 

2kik2k3 fkkik2 1 
= niQ ^ h 



aq^ \ 2q^ 

The difference of these two equations gives: 
'kkik2 



A=l 



0. 



(4.21) 



p+- p- 



2g2 



+ J 



1 1 



+ 



0. 



(4.22) 



This is precisely equation (I4.19p and this observation is the analog of the fact that, when the 
charge densities are constant, the sum of the bubble equations vanishes. One can write (I4.22p as 



2q(Y,kAQ^i ( 

^ A=l ^ ^ 



kkik2 
2g2 



+ J . 



(4.23) 



This equation determines the scale- invariant ratio, r'^jr'_ while the remaining bubble equation 
in (I4.2ip determines vhq. 

Thus, once one specifies all the charges, the geometry is fixed. The non-trivial constraints 
come from the fact that the charge densities on the supertube are all related and thus the charges 
Qa and J must be related. Indeed, (I4.16P gives: 



4:7rqp2{ip) 



kki 



, , , kk2 C: 

47rgpi(V') + —F 



-1 



A'Kqpi{il))F + 



kko 



(4.24) 



One then gets p(^/') from (]4.17p . In particular, from (I4.18P one has 



J 



ki f^^'i 
Tq 



/ , N kk2 C; 

4:Tiqpi{'4)) + F 



4:7Tqpi{ip)F + 



kko 



dip 



(4.25) 
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5 Comparison with probe supertubes 



In the previous section we obtained the regularity conditions near the wigghng supertube in the 
fully back-reacted solution with two GH points and a supertube. Obviously one would like to 
generalize this solution to include more supertubes and more GH points but, as we discussed 
earlier, such fully back-reacted solutions are, at present, beyond the limit of explicit computation 
and so we must seek another approach to this problem. 

It is natural to fall back upon a probe approximation and use the DBI action that describes 
supertubes in the regime of parameters in which they do not back-react on the geometry. This 
action can be used to analyze more general configurations, but the problem then becomes one 
of determining the extent to which the DBI action captures the correct physical behavior of 
the complete physical (back-reacted) solution. In [19], it was shown that for supertubes with 
constant charge densities, the DBI action of supertubes always precisely captures the conditions 
necessary for smoothness and regularity of the fully back-reacted metric at the supertube location. 
Our purpose here is to generalize this to supertubes that have non-trivial density modes, and 
show that the solutions of the DBI action always correspond to supertubes whose back-reacted 
solution is smooth and free of CTC's. We show that not only can one use the DBI action to 
extract the local regularity around the supertube but one can combine the information coming 
from this action with the general analysis of solutions in section H] to infer the full regularity 
conditions at the GH points that lie far from the supertube, and thus find the complete set of 
functional bubble equations that govern multi-center solutions containing one wiggly supertube 
in an arbitrary Gibbons-Hawking base space. 



5.1 Regularity conditions from the supertube Born-Infeld action 

The Born-Infeld analysis of supertubes is best done in the duality frame where the supertube has 
DO and Fl electric charges and D2 magnetic dipole charges, and for this we first have to dualize 
the background given in (12. ip and (12. 3 p to this type HA frame. The details of this transformation 
may be found in [19] and here we will also adopt the same notation and conventions. 

We first reduce to ten dimensions along xs and then do two T-dualities along x-j and xg. The 
Fl string and the dipole D2 brane are wrapped along the xq = z. The resulting metric (in the 
string frame) is: 



dsi, = -^=l=-idt + kf + ^/Z^zldsl + ^^^dz' + \^dsl, (5.1) 
and the dilaton and the Kalb-Ramond field are given by: 

e* = — A_ ] B = (Zf ^ - l)dt Adz + Z^^k A dz . (5.2) 

\Z3Z1 J 

The non-trivial RR potentials are: 

= {Z^^ - l)dt + Z:^^k , (5.3) 
C^''^ = [Ca + V-^K^if^) n'L^ Adz- {Z^\dt + k)A 5(2) + dtA A dz , (5.4) 
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where k and Zj are given respectively by (EH]), f l2Jl]) . fl2A4ll and f l212ll . and 



V X C 



-VL, . 



(5.5) 



We will denote the world-volume coordinates on the supertube by and = 6. To 

make the supertube wrap z = Xq we take C,^ = z and we will fix a gauge in which = t. Note 
that z G (0,27rL2). The profile of the tube, parameterized by 6, lies in the four-dimensional 
non-compact Gibbons-Hawking space. We will here take the profile to be round, 6 = kip, with 
dipole charge k, but will allow for density modes: the gauge field, J-", living on the world-volume 
of the supertube will depend on 9. We take 



2TTa'F = T = J^tzdt Adz + J^zedz A d6 . 



(5.6) 



Supersymmetry requires that J'tz = 1 [18], but J-'ze, which corresponds to the DO-charge density, 
can be an arbitrary function, of 9, Tzq = J^ze{9). 

The supertube action is a sum of the DBI and Wess-Zumino (WZ) actions: 



S = -Td2 j d^^e-^^-det (Cab + B^b + J'ab) + Td2 j d:'i[C^^^ + C^^^ A (^ + B)] 



(5.7) 



where Gab and Bab are the induced metric and Kalb-Ramond field on the world-volume of the 
tube. We have also chosen the orientation such that e^g = 1, and throughout this section all the 
functions will be evaluated at the location of the supertube, and we will omit the hats that we 
used to note this in section H] 

After some algebra, the action simplifies to: 



S = Tn2 I d'^ 
1 



1 \ 



nl/2 



[{K' - V(^(l - Ttz) - J'ze)? + VZ^Z^il - J-t,)(2 - Zi(l - J'tz))] 



For a supersymmetric configuration {J^tz = 1) we have: 



. (5.8) 



Sj^t^=i — Sdbi + Swz — ~Td2 / dtdzd9 J-' ze ■ 



The foregoing supertube carries Fl and DO "electric" charges, given by: 



ST 



Fl 



d9 



dC 



tz 



N. 



ST 



D2 



DO 



dzd9 Tzd ■ 



(5.9) 



(5.10) 



To avoid the proliferation of parameters, it is convenient to use the system of units described in 
the appendix of [T^. With these conventions, the supertube supergravity charges, Q\ and Qi-, 
and dipole charge, /c, are related to the integer charges, iVf^, N^^ and k via: 



-AT' 



ST 

A 5 



k = -qk, 



(5.11) 
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and the values of the tensions are: 

Too = 1 , SttL^ Tpi = 1 
The Hamiltonian density is: 



27rTi 



D2 



Tl 



Fl 



tz 



(5.12) 



(5.13) 



where we introduced pi and p2, the charge densities of Fl and DO charges respective!}!^: 



Pi 
k 



1 dC 



tz 



21 
k 



D2 



dz 



zQ 



(5.14) 



Note that, in the Born-Infeld calculation, we are treating a supertube of dipole charge A; as a 
supertube of dipole charge one wrapped k times, and hence the ^^ are densities per unit strand 
of the supertube, and the total density is pi. We want to remind the reader that pi, p2 are not 
constant but vary with ^l). 

One can easily integrate this to get the total Hamiltonian of the supertube 



dQ nw,^=i 



ST 



(5.15) 



Thus the energy of the supertube is the sum of its conserved charges which shows that the 
supertube is indeed a BPS object. 

In the DEI analysis, the first equation determining the regularity of the supertube comes 
from the relation that gives the Fl charge density pi 



Pi 



dC 



tz 



kTi 



D2 



4T, 



Fl 



+ 



ze 



(5.16) 



Using the relation (15. lip for /c, the values of the tensions fl5.12p and the definition of p2 (I5.14p . 
this can be rewritten as 



47rg pi + k- 



v 



Anq P2 + k 



V 



k 



2 ^3 

V 



(5.17) 



and this is exactly the same equation as ( 14.16^ . 

The second regularity condition comes from expressing the angular momentum of the tube 
as a function of the charge densities. Note that, in contrast to the supergravity construction in 
which this angular momentum is a free parameter, in the DBI calculation the supertube angular 
momentum is related directly to the other charges. The angular momentum density p along the 
ip circle is, in general, given by 



k 



dip 



(5.18) 



^The factor of 4 in (|5.13p comes from the relation (|5.1ip for Q 



A- 
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One should note that, even if one is interested in time- independent solutions, in order to compute 
p one has to allow a time dependence initially, compute and p, and only then impose time 
independence. The angular momentum of the supertube is then given by 



J = ddp. (5.19) 



After lengthy, but straightforward computation, one finds: 



P 



' ( - ^^K^ - kV-^^^^ I , (5.20) 



Stt^ \^ k 4TTqp2 + k^ 

where, in the second line, we have used the relations (15. lip . (I5.12p and (I5.14p . along with the 
fact that J^ze is independent of z. Using (I5.17p . (I5.20p can be rewritten as: 

47rq{2pV + piK^ + pa^') = kLs , (5.21) 

or again 

kp = 2Tiqpip2. (5.22) 

Since the densities pi, p2 and p are functions of ip, these two equations are functionnal equations, 
depending on tp. They exactly match (14.151) and (14.171) . Hence the supergravity and the Born- 
Infeld calculation give the same regularity conditions for supertubes of arbitrary charge density. 



5.2 Obtaining the bubble equations from the DBI action 

In the last subsection, we showed that one can obtain the regularity conditions at the supertube 
location directly from the DBI action, even for a wiggling supertube. Since the supertube is 
treated as a probe, the regularity conditions for the fully back-reacted solution at the other GH 
points - in other words the bubble equations - are not captured by the DBI analysis. Nevertheless, 
we will show here that the structure of the general solution outlined in Section [2], and in particular 
the regularity conditions at the GH points allows us to infer the complete set of functional bubble 
equations starting from the regularity condition at the supertube location. 

In order to obtain the bubble equations from the regularity conditions of a supertube with 
variable charge densities, we need to first understand these equations for a supertube with con- 
stant charge density in a multi-center GH solution. Since this solution does not depend on the 
Gibbons-Hawking fiber, it descends to a multi-center four- dimensional solution of the type first 
constructed in p2] , where the GH centers are D6 and anti-D6 fluxed branes and the supertube is 
a fluxed D4 brane. For a solution with GH centers and one supertube, there are bubble or 
integrability equations insuring absence of Dirac-Misner strings at the GH points, and one such 
equation for the supertube point. The sum of the (A^ -|- 1) bubble equations is zero. 

An important feature of these equations is that they only depend on the pairwise interaction 
between the centers, and do not contain terms that depend on the product of charges from 
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three or more points. For example, each term in the supertube bubble equation depends on the 
supertube charges and on the charges at one of the GH centers. There is no term proportional 
to the product of two or more far-away GH points. 

As we have shown in the explicit construction of the fully back-reacted solution, the two 
properties: the fact that the bubble equations sum up to zero and only involve pairwise inter- 
actions, are still true for solutions with wiggling supertubes. It is also easy to see that these 
two properties hold more generally for solutions with arbitrary number of centers and varying 
charge densities. Recall that the bubble equation at a given center can be obtained from the 
requirement that there be no Dirac-Misner (DM) string starting at that center. Suppose one has 
already ensured the cancellation of DM strings at all points but one. Then the last point can 
only have a DM string ending at infinity and if one has already imposed regularity at infinity, it 
follows that the last bubble equation is automatically satisfied. This establishes the first property 
of the bubble equations, that is, their sum is zero. 

The pairwise property of the bubble equations follows from the general structure of equation 
( I2.23p . determining the angular momentum vector u: at both the GH points and at the supertube 
point, the right-hand side of equation fl2.23p is sum of products of the local charges with harmonic 
functions that satisfy linear equations. Hence, the contributions from various centers to the 
equation for the supertube add linearly. 

In order to write the bubble equation for the supertube in the usual form, one should observe 
that equations f l5.2ip and f l5.22p imply that if the supertube is round and hence it sits at a point 
in the base of the GH space, the two electric densities and the angular momentum density 
are not independent. Hence, this supertube has one function worth of degrees of freedom. 

If one now integrates equation fl5.2ip over the supertube world-volume, one obtains the stan- 
dard bubble equation for the supertube: 

QiKilsT + Q2K2\sT - kL^lsT + JV\sT = , (5.23) 

where Qi and Q2 are the total electric charges, and J is the total angular momentum of the 
wiggly supertube f l5.19p . 

Hence, one can imagine taking a round supertube of given total charges Qi and Q2 and letting 
the charge densities wiggle in accordance to fl5.2ip - fl5.22p . The first three terms in the supertube 
bubble equation f l5.23p remain the same, and only the term proportional to the supertube angular 
momentum, J, is modified. Since one assumes the total charges to be fixed, it is clear that because 
there is only one independent density function, the location of the tube will change in order to 
satisfy (ESSD- 

If one now remembers that the Ki, K2, and V harmonic functions are sums of harmonic 
functions sourced at each of the Gibbons-Hawking points, one can easily see that the supertube 
bubble equation is entirely composed of terms of the form: 

Qik{ 
\rsT - n 

plus, possibly, some constants. 

Using the fact that all the terms in the bubble equations for the other GH center contain 
pairwise interaction terms, and the fact that the bubble equations sum up to zero, we find that 
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(5.24) 



the only effect of the wiggly supertube in the bubble equations for a given GH center, is to 
modify the terms 



Jq' + Qik{ + Q2k\ - kQ\ 
I^ST - ri\ 



(5.25) 



by replacing J, Qi, Q2 and k with the corresponding quantities for the wiggly supertube. Further- 
more, given that k does not oscillate, and that one is considering a supertube of constant total 
electric charges, the only modifications brought about by the wiggles show up in the expression 
of J. 

Hence, the functional bubble equations for backgrounds containing a wiggly supertube in 
multi-center GH solutions are exactly the same as those for a round supertube, except that one 
replaces the angular momentum J of the round supertube with the angular momentum of the 
wiggly one, which is obtained from the density modes via equations fl5.19p . fl5.22p and f l5.2ip . 



6 The entropy enhancement mechanism 

Having constructed the full solution corresponding to a wiggly back-reacted supertube in an 
ambipolar space, we now turn to the exploration of the amount of entropy such a supertube 
can store. In the regime of parameters where one can ignore the back-reaction of supertubes, 
this question was studied using the Born-Infeld action [25], and it was found that the entropy a 
supertube can store depends on its location in the solution, and is much much larger than that 
of a supertube in flat space. In fact the entropy yielded by the Born-Infeld action diverges as the 
supertube approaches the hypersurface between the regions of opposite signature in the ambipolar 
GH base. It is therefore crucial to understand what are the limits of entropy enhancement, what 
cuts it off in the fully back-reacted solution, and what is the dependence on the total charges of 
the entropy that a supertube can store. 



6.1 A near-critical supertube 

The maximum entropy enhancement of a supertube potentially occurs when the supertube ap- 
proaches the critical surface, = 0, or = r_. In this limit, \F\ (given by (14. 8p ) is large. 
For generic values of the charges (such that Y^a=i ^aQa 7^ 0), equation (14.231) implies that the 
supertube angular momentum J becomes large in this limit, of order \F\. Moreover, the large 
|F| limit of 04.251) gives 

J = jT- — / piiip) (6.1) 

k K2 ■''0 

Hence, in r+ — r_ limit both J and the density fluctuation, J^^'^ PiiS^Y dip, grow proportionally 
to \F\. In this limit, equation (I4.24p becomes 

/ n ^1 , WfclP?(^) 

k2 kkl F 
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note that the term proportional to remains finite in the hmit, and hence it cannot be 

discarded. 

We would like to argue that (16. ip means that J must be negative, and since we have taken the 
ki to be positive, this means that we need to argue that k is necessarily positive. The simplest 
way to see this is to consider the scale factor of the Kaluza-Klein circle in the D1-D5-P duality 
frame: 

^3 K^K2 + L^V 



W 



y^Z^2 ^/{K2K3 + LiV){KiK3 + L2V) ' 



On the critical surface, one has 



W 



k.+ 



k r± 



2A 



(6.3) 



(6.4) 



If k and ^3 have opposite signs then there is a danger that the equation ^ = — ^ will admit 
a solution and when this happens the warp factor W diverges and the metric becomes singular. 
We will therefore require that all the kj and k are positive. 

This means that, in the near critical limit, J is necessarily negative and proportional to the 
mean-square of the fiuctuation density. If the supertube has large density fiuctuations then J is 
necessarily large and negative. 

6.2 A simple example 

We now illustrate the results above by computing a very simple example. Take pi to be given 
by: 

Qi 



47rg 



— (1 + a COS'?/') 



(6.5) 



that is, we only take the first Fourier mode. The other densities, p2 and p, are then fixed by 
regularity as in (I4.16P and (I4.17p . One can integrate (I4.24p to get 



- kki ^ k'^kik2 
Q2 = -r + 



q 



{F^ - 1) 



(6.6) 



which may be rewritten 



4 + ^^ii? 

q 



k?'k\k2 

g2 



(F- - 1) 



(6.7) 



One can then get J from f l4.22p 

J = 



kkik2 
2g2 



(6.8) 



For fixed Qa, (16.71) means that F, and hence, r'^/r'_ must be readjusted as the amplitude 
of the oscillation grows. For the sake of definiteness, assume that the Qa as well as the kj and 
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q, k are all positive. In order to have a large entropy, we need to get a large, negative angular 
momentum. This is done by taking r_|_ — r_ — > 0_ so that F is very large and positive. In this 
limit, one can see that has to scale like F in order to satisfy (16. 6p . One can thus assume in 
this limit, without loss of generality, that 1 « a « F. Combining then (16. 7p and (16. 8p one 
finds: 



J 



Q1Q2 
2k 

Qi Q2 
2k 



a' 



a 



01 

Ak 



Qi + 



ki 
Akk, 



kko 



-1 



Q2 



kki 



+ 



\ p2 



(6.9) 



Thus taking J very large and negative corresponds to extremely large fluctuations, a, and these 
can be sustained by taking the supertube extremely close to the critical surface (a^ ~ F). 



6.3 Large-distance behaviour and an angular momentum bound 

If there is a large angular momentum then there is necessarily a limit on this given by the 
requirement that there are no CTC's at infinity. Since the supertube wraps the if) circles, we 
would expect the primarily limitation would come from this direction, and for fixed t the metric 
coefficient of {dip + is: 

Y = - {Z.Z^Z.y^/' + {Z,Z2Zsy/'V-\ (6.10) 

and this must be non-negative everywhere. 
At infinity (r — )■ 00) one has: 

Y "I^^ r, (6.11) 

(kMh + \k)Y" 

from which it follows that mg must be negative. It then follows from the bubble equations (14.211) 
that 

kikikz fr'+ , r'_\ kifcs o 1 ^' , ^ n (a^o\ 

+ J = -imQ[r^ + r^) > 0. (6.12) 



\ a a / 2q'^ 

This places a bound on how negative J can be. Indeed, if the supertube approaches the critical 
surface, then J, and hence the mean-square fluctuations, are limited, not by the intrinsic charges 
of the supertube, but by the product, kik2k2, and thus by the total charge of the corresponding 
black hole or black ring. 

One can also use (I4.23P to write (I6.12p as a bound on how close the supertube can come to 
the critical surface: 

n \ ^ ^ 2kik2k3 
^ / (r_ - r+) q 

Note that for the supertube to come extremely close to the critical surface the fluctuations 
must be large so that J is large and negative, and so (I4.23P means that F should be large and 
positive, where one should remember that to arrive at this we have assumed that r'_ — r', > 0. 
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It is important to stress that the DBI analysis, that only captures the local properties of the 
back-reacted solution, cannot encode large distances behavior, and hence the bound on J. The 
knowledge of the fully back-reacted solution was therefore crucial to obtain this bound. 



6.4 How to get most entropy from a supertube 

Recall that the entropy of the supertube is given by 



S ^ ^QiQ2 - 2kJ. (6.14) 

In flat space A; J is positive and proportional to the radius of the supertube and thus the angular 
momentum is bounded according to |J| < \^QiQ2\, where the Qa are the intrinsic charges of 
the supertube. On the other hand, it was shown in |25l[19], using a brane-probe approximation, 
that in deep scaling solutions, the dipole-dipole interactions could allow kJ to become arbitrarily 
negative and therefore a supertube could store a vast amount of entropy, far beyond its flat-space 
limit. This is the entropy enhancement mechanism. 

We have now seen precisely the same process in the fully back-reacted solution, and how 
fluctuations in the deep-scaling limit of an AdS throat can indeed lead to very large negative 
values of kJ through dipole-dipole interactions as the supertube approaches the critical (V = 0) 
surface. The advantage of the back-reacted solution is that we can also see that there is a bound, 
(I6.12p . on just how negative J can become. If F is large, one has ^ r'_ and the supertube 
electric charges become irrelevant. Thus the bound becomes: 



kik 



- J < ^{^lk3 + k). (6.15) 



2q 

r' 

where 7 = — defines the aspect-ratio of the triangle defined by the supertube and the GH 
points. If the flux parameters, kj, are very large, then one has the following bound on the 
entropy enhancement: 



S ~ \/QiQ2 - 2kJ < ^lQ,Q,+ '^{4^k, + k) ^ ^^(47A;3 + fc). (6.16) 

Thus the entropy enhancement can be extremely large but is still limited. In particular, the 
limit involves precisely the dipole-dipole interaction between the supertube and the background 
geometry. In fiat space the supertube entropy is limited by its intrinsic charges, whereas here it 
is limited by the charges of the complete background. We have thus demonstrated that entropy 
enhancement is a very real, and potentially very large phenomenon but that it is bounded. 

For one supertube, one can use the fact that the M2 charges of our solution are proportional 
to the product of the magnetic fluxes on the cycle between two GH centers to estimate the 
dependence of the entropy on the total charges. Assuming that the supertube dipole charge k is 
of the same order as the fcj, this gives S ~ \/k^ ~ \fQ^ . This entropy has the same growth with 
charges as that of a normal supertube in fiat space; the difference is that it is not realized by 
putting all the charges on a supertube, but rather by using the charges to create a two-centered 
bubbled solution, and putting a very small supertube in this background. 
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Hence, the entropy of a single supertube in a bubbled background, though very enhanced, 
does not give a parametrically-larger entropy than that of a two-charge system. To get more 
entropy, we have to find a way to put a larger amount of negative angular momentum on this 
supertube, without destroying the asymptotics of the solution. It is not hard to see that the 
angular momentum bound (16.151) is much like that of a BMPV black hole: the M2 charges of 
the background are proportional to the product of the fluxes: Qi ~ ^2^3, and the bound of the 
absolute value of J is 

\J\ < VQ1Q2Q3 , (6.17) 

exactly like the BMPV black hole. However, it is well-known that for a given set of charges, 
the BMPV black hole is not the object with the largest angular momentum. A black ring, or a 
supertube can carry a parametrically-larger angular momentum, proportional to the square of 
the charges [IHl SS] • 

Thus, a possible way to obtain more entropy than \/Q^ is to place more than one supertube 
in the two-center bubbling solution. Indeed one could use a large non-wiggly supertube to act as 
an angular momentum sink and another smaller one that will give the entropy via the entropy 
enhancement mechanism. The two supertubes could be given dipole charges that are oriented in 
the same direction so that there is no danger of W (equations fl6.3l) and (16.41) ) vanishing. However, 
the angular momentum of the non-wiggly supertube will point in the opposite direction from that 
of the entropy-enhanced supertube. Assuming that the dipole charge of the wiggly supertube is 
again of order Q^^'^, its angular momentum is now of order Q^, and hence the total entropy of 
the system is 

S \lk\J\ ^ (6.18) 

while this entropy is still not black-hole-like, it is parametrically larger than that of the two- 
charge system. Of course to establish that this kind of entropy is indeed present in our system, 
one would need to construct the fully-back-reacted solution with a wiggly and a non-wiggly 
supertube in global AdS^ x 5*^ and to check that this metric is free of CTC's. 

Another possible way to obtain a larger entropy is to place many wiggly supertubes in this 
background. In general, the dipole charges of these supertubes should have the same orientation 
as the ^3 of the background, in order to avoid zeros in the scale factor W in equation (16. 4p . 
Hence, when they spin, their angular momenta will be oriented along the same direction, and 
will sum to the total allowable value of \ J\ (which is of the order Q^^'^). The total entropy will 
contain, besides the entropy of each individual supertube, a component coming from the many 
ways of partitioning J between various supertubes, and will also be generically larger than \/Q^. 
Of course, by adding a large non-wiggly supertube that acts like an angular momentum sink we 
can increase the upper limit of | J|, and presumably obtain an entropy larger than even (I6.18p . 

Probably the best way to bypass the angular momentum bound and obtain a large entropy is 
to put two counter-rotating wiggly supertubes, which naively could have huge opposite angular 
momenta, and could give us as much entropy as we want. However, there is no "free lunch:" 
the dipole moments of the two supertubes now have to be opposite, and for most of the possible 
locations of the supertubes, this will cause problems for the scale factor W in equation (16. 4p . 
However, it may be possible to select relative locations for the two supertubes in the vicinity 
of the critical surface such that W will never have a zero. If such a configuration exists, and 
satisfies the bubble equations, one can imagine starting to wiggle the supertubes, in order to 
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get more and more states. As one does this, the bubble equations determining the relative 
positions of the points change, and one can imagine that at a certain upper value of the angular 
momenta these positions will become incompatible with an everywhere-positive W. It would be 
certainly interesting to explore this configuration in detail, and to see how much entropy can two 
oppositely-spinning supertubes store, and how does this entropy compare to that of a black hole 
with the same charges. 

7 Conclusions 

We have constructed an infinite-parameter family of smooth supergravity solutions that have 
three charges and three dipole charges, and that are microstates of black holes of classically- 
large horizon area. These solutions are obtained by placing a two-charge supertube of arbitrary 
shape in a two-center Gibbons-Hawking base space. In general the solutions for such supertubes 
can be implicitly written using scalar and vector Green functions, but so far no explicit solution 
has been written down, owing to the complicated form of the Green functions, and to the fact 
that for most of the base spaces that are physically of interest, the ambipolar GH spaces, the 
Green functions were unknown. 

In general, the oscillations of a supertube in the four-dimensional base space of the solution 
can be parameterized by four continuous functions but we have identified a subclass of super- 
tube oscillations where the shape of the supertube remains round, and only the distribution 
of electric charges inside the supertube world-volume changes. The solutions corresponding to 
these supertubes depend on one arbitrary continuous function, and their magnetic dipole fields 
are exactly the same as those of round supertubes, which considerably simplifies their explicit 
construction. The main ingredient that enters in the construction of the explicit solution is the 
scalar Green function on the base space. Given that the ambipolar GH base spaces have regions 
of signature —4 and +4, with intervening "critical surfaces," the scalar Green functions is much 
more complicated than it is for regular GH spaces. We could only find its explicit form for the 
two-center ambipolar Gibbons- Hawking space given by the harmonic function j^^^ — j^r^- This 
was done via a highly non-trivial procedure that involved reducing the five- dimensional Green 
function on the smooth Lorentzian space-time that can be constructed from this base, that is, 
from the Green function on global AdS^ x S"^. 

In constructing this family of smooth horizonless black-hole microstate geometries, we have 
also found that the bubble (or intcgrability) equations that determine the relative locations of the 
Gibbons-Hawking centers and of the supertube are unchanged for the GH centers but become 
non-trivial functional bubble equations on the supertube. Remarkably, the same functional 
bubble equations can be recovered by examining the Born-Infeld action of a probe supertube 
in this space. Given the completely different nature of the two calculations that yield the same 
functional bubble equations, and given that one calculation is done in the regime of parameters 
where the supertube does not back-react on the geometry, while the other was done in the regime 
where it does, this result has several important implications. First, it points to the existence of 
a non-renormalization theorem that protects the functional bubble equations as one moves in 
moduli space. Such a non-renormalization theorem exists for the four-dimensional multi-center 
solutions that come from a five-dimensional solution with a tri-holomorphic U{1) invariance 
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(and hence a Gibbons- Hawking base space), and our analysis finds that this non-renormahzation 
extends to solutions that do not have this invariance. 

One of the most important uses of the non-renormalization theorem for these [/(l)-invariant 
solutions has been the quantization and counting of certain finite-dimensional moduli spaces 
of multi-center configurations [HI E]. This is done at weak coupling, and then extrapolated 
to the regime of parameters where all the branes back-react using this non-renormalization 
theorem. The fact that the infinite-dimensional moduli space of wiggly supertubes does not 
receive quantum corrections implies that if one quantizes this moduli space at weak coupling 
(using for example the Born-Infeld action of supertubes) one can find how much entropy comes 
from fully-back-reacted wiggly supergravity solutions. 

Another important result of the agreement between the DBI analysis and the fully back- 
reacted descriptions of the bubbling solutions containing supertubes is that one can streamline 
the construction and analysis of supertube configurations that give smooth microstates without 
constructing the full supergravity solution for each and every microstate. Indeed, as shown in 
Section m and in [Tn|, the local condition that the supergravity solution near the supertube be 
smooth and free of closed timelike curves is exactly the same as the condition that the supertube 
be a solution to the DBI action. In addition, we have also seen in this paper that the conditions 
that the solution be free of closed timelike curves near the GH centers, that is, all the other 
bubble equations, can also be obtained from the functional bubble equations for the supertube. 
Hence, if one wants to find the properties of a solution containing a supertube that has a charge 
density given by an arbitrary function pi{6), one first determines the density distribution for 
the other charge P2{0), using equation (14.161) . then finds the total supertube angular momentum 
density J (9) using equation (I4.17p . The next step is to integrate these densities to find the total 
charge and angular momentum, and to use the bubble equations to determine the location of 
this supertube and of the other Gibbons-Hawking centers in the full solution. Our analysis finds 
that as far as the full supergravity solution is concerned, the smooth supertube with variable 
charge and angular momentum density behaves exactly as a singular f/(l)-invariant supertube 
that has the same total charges and angular momentum. Hence, to analyze whether a given 
wiggly supertube gives a smooth and regular geometry upon back-reacting one simply has to 
construct the fully back-reacted solution of an equivalent round supertube, which respects the 
tri-holomorphic U{1) invariance of the GH base, and which is straightforward to write down in 
terms of harmonic functions. 

Last, but not least, our analysis finds that the entropy enhancement mechanism, first un- 
covered in [25] for non-back-reacted supertubes, extends to fully back-reacted solutions. This 
mechanism allows supertubes with relatively small electric charges to have a much larger entropy 
in a background with large magnetic fields than in flat space. We showed in Section E] that the 
entropy of such a supertube is bounded above not by its own electric charges, but by the electric 
charges of the background in which it is placed. This establishes that the entropy enhancement 
mechanism is not an artifact of the Born-Infeld approximation to the supertube dynamics, but 
is a feature of fully-back-reacted solutions containing supertubes. 

Clearly, the most important open problem raised by our work is to determine how much 
entropy can be found in the class of solutions that we have constructed. In contrast to all three- 
charge solutions that have been obtained so far, whose moduli space is finite-dimensional, the 
solutions we construct have a much larger moduli space, whose dimension is infinite. Hence 
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they should have a much larger entropy than that found by the semi-classical quantization of 
the highly symmetric bubbled solutions in fU] [Ti] . The question is how much more can be 
gained through the density fluctuations and through entropy enhancement. Indeed, the entropy 
of [/(l)-invariant (and toroidally invariant) microstates solutions is parametrically smaller than 
that of the corresponding black holes, as one might have in hind-sight expected from counting 
only microstates that respect a certain isometry. It is also clear that the solutions we construct 
are not the most general black hole microstates one can flnds in supergravity. All of our solutions 
are independent of the internal space (which one can take to be or a more general Calabi- 
Yau) , and utilize only one of the four functions worth of oscillations that a supertube can have in 
Nevertheless, from the Born-Infeld action of supertubes we know that their moduli 
space is characterized by continuous functions, and the isometry we break does not destroy the 
nature of the moduli space and the nature of the counting problem, it just uses one of the 
eight bosonic degrees of freedom in the DBI action. Consequently, we expect the entropy of 
the solutions we construct to differ from the entropy of the supergravity solutions coming from 
all the possible supertube oscillations by a numerical factor of \/8 (or \/l2 if one considers the 
fermion partners of the bosonic oscillations). Hence, counting the solutions we have constructed, 
while not capturing all the entropy of supergravity solutions, will give a flnite and known fraction 
of the complete set of solutions. In contrast, the extra U{1) isometry of the solutions counted 
in [m [H] destroys the function-dependent nature of the moduli space and makes it flnite- 
dimensional, hence such solutions cannot be expected to have an entropy that grows with the 
charges in the same way as that of the most general smooth supergravity solution. 

On the more technical side, one drawback of our work is that we have not been able to obtain 
the expression of the four- dimensional rotation parameter u in closed form. The only danger 
associated with co is that it could lead to the appearance of closed timelike curves, either via 
Dirac-Misner strings, or more globally. The absence of the former is guaranteed by the functional 
bubble equations. Furthermore, as we have argued above, all the global properties of a solution 
with smooth wiggly supertube can be captured by a Gibbons- Hawking solution that contains an 
equivalent round supertube. If the co of the latter solution does not cause closed timelike curves 
far-away from the supertube location, it is likely that the u of the wiggly solution will not cause 
problems either. This being said, it would still be quite interesting to try to obtain u in closed 
form, at least for some of the solutions, and see that indeed the regularity of the solutions is 
insured by the functional bubble equation and by the regularity of the solution with an equivalent 
round supertube. Finding u explicitly would also probably be essential to the complete analysis 
of a solution with multiple counter-rotating, wiggling supertubes. We are currently trying to 
develop a class of solutions that contain such supertubes. 

The work presented here has demonstrated the reality and viability of the entropy enhance- 
ment mechanism and thus represents signiflcant progress in flnding microstate geometries. In 
subsequent work we hope to be able to construct solutions with multiple wiggling supertubes and 
show how they can interact with one another and generate mutual entropy enhancement while 
keeping the total angular momentum small. The interactions and closed-timelike-curve analy- 
sis will limit the enhancement but we remain hopeful that it will get us near the long-sought, 
semi-classical black-hole entropy. 

^For two-charge solutions that utilize the internal and fermionic degrees of freedom see [T71 H7] . 
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Appendix A. Details of the Green function 

In the main body of the paper we only considered the details of the Green function, f l3.50p . 
that were directly pertinent to finding the wiggling supertube solution. In this Appendix we 
will examine the Green function in more detail and show that it does indeed have the requisite 
properties. 

Appendix Al. Branch cuts and singularities 

We begin by considering the possible singularities of G and while the physical domain has \z\ = 1, 
we will consider \z\ < 1 for the same reasons that were outlined in Section [3.61 We will show 
that the only singularity is precisely the physical one exhibited in Section 13.7.11 

Prom (I3.45P it follows that y/a2 + ib2 vanishes only if \u — vz\'^ = \yz\'^ while from fl3.44p one 
has \u — vz\ > 1, \yz\ < 1 with equality if and only if z = 1, ^ = ^ = = 7i/2. Using this 
one then sees that 02 = 62 = if and only if ^ = C,', ip = = 0' and 6 = 6' = it /2: In other 
words both points coincide with each other and lie on the critical {V = 0) surface. Thus the 
only potential singularity generated by y/a2 + 262 occurs only for a particular class of coincident 
points and this singularity was analyzed in 13.7.11 and we will also discuss some further details 
below. More generally, from the fact that \u — vz\ > 1, \yz\ < 1 it is not hard to convince oneself 
that 02 is everywhere non-negative and so there is a globally analytic, single determination of 
the square root, which means that there are no branch points and branch cuts in the region of 
interest. 

Now consider the singularities associated with the vanishing of the function, A. From the last 
expression in f l3.59p we see that there are possible singularities in when either a) A = 0, b) 
U = T = or c) 6 = 7T — 6' . Note that while |f/| — )■ 00 is certainly possible if r± + r'_^ = A, the 
last expression in fl3.59p shows that remains finite. The other expression in (13.590 . along with 
(I3.57p . shows that there is, in fact no singularity associated with A vanishing by itself. Similarly, 
the same expression, combined with the form of m in f l3.46p . shows that there is no singularity 
associated with 6 = 71 — 6' alone. Thus the only singularity in fl3.59p occurs when U = T = 0, 
indeed there is a double pole if one sets T = first and then takes U 0. We will now show 
that this double pole is generically cancelled by the other terms in the complete Green function 
and that the only singularity arises when the points are coincident. 

Setting T = in the Green function is equivalent to taking z = f^. In terms of the earlier 
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notation, this means that: 

^ ^ (^+ - r) 

where we have used the fact that z+z^ = ^. The possible double pole thus emerges in the limit 
z+ ^ Note the order of limits here: The value of z is fixed first and then the limit of z+ is 

\m\ 

taken. Also note that the choice z = has |z| = 1 and is thus in the physical domain. This 
means there is the possibility of strongly singular behavior at non-coincident points in G and we 
will now show that such singularities do not occur. 

We begin by assuming A 7^ and hence V 0, then the limit of interest can only happen 
if + x)"^ — y"^ = {cos6 + cos 6*') = 0. Thus the potential singularity can only occur on an 
"image surface" at the "conjugate latitudes," = tt — on the S^. 

Let 6 = TT — 6' ~ e and so + x)^ — ~ ^{^)- O^is can also easily verify that 

^ = V-\ = + ' = — ± —sineVV + Oe^ , A.2 

\m\ ue + V \m\ 2m 

and therefore {z+ — j^)^ ~ 0{e^). On the other hand, for physical parameter values = 1) 
and non- coincident points, one has 02 + ^62 7^ but one also has 

u = - (i + x){{u + ve-''') - z{ue-''' + v)) + {{1 + x) - y){v - uz)e-''' ~ ^(e^) , (A.3) 

where we have used (]A.2p and ((1 + x) — y) = cos^ ^{9 + 9') ~ 0{e'^). It follows that u vanishes 
strongly enough to cancel the potential singularity as e — )■ and that the Green function fl3.50p 
is, in fact, finite as e — i- 0. 

Now suppose A = 0, then (,=^',9 = 9' and = 0', then m = ^-^^ and z± = 1. One then 
finds that: ^ ^ 

A " m{z-iy ' ^^'^^ 

which indeed has a double pole as 2; = e^^^~^ ^ — > 1. This is, of course, the correct double pole 
for the Green function as ip ^ ip' with = C,', 9 = 9' and = 0'. 



Appendix A2. The scalar fields sourced by simple sources 

In this appendix we compute the moments of the Green function by integrating it against the 
Fourier modes of a given electrical charge distribution. By putting together these moments one 
can reconstruct the full warp factor sourced by supertubes with variable charge distributions in 
the two-center ambipolar space discussed in this paper. 

Fourier modes along the fiber 

Taking g = 1 for the present, we compute the moments of the propagators, that is, the integrals: 

= ^di^ e^^'^'^') G = -l-z^G. (A.5) 
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where the integral is over a circle of radius 1. We write the propagator as 

^ 0, A + 102 

We assume m > 0. The integrals with m < can be derived from the reality condition 

I-m = rm- (A.7) 

The function (j){z) is regular for \z\ < 1 but has a branch cut in the region \z\ > 1. Then, to 
compute the part of the integral containing it is convenient to change the coordinate 

z -)■ z~^. 

I^ = -l^z^ <j>{z- x) + -(f^z-^ ^{z- -x) (A.8) 
I J aJ\ I J aA 

where 

A = z^A{z~^) = m{z - z-^){z - zZ^) • (A.9) 

The first integral in ( 1A.8I) receives contributions only from the pole of A inside \z\ < 1, which is 
z_ if ?7 > and z+ if U > 0: 



^ (p ^z"" (t){z- x) = 27r _ , , if U^O. (A.IO) 



dz ^ ^. , ^ cos ^ + cos 6' 

— z^ (j){z; x) = 27r 

a A am[z+ — z^) 

Using the identities derived above, one can re-write this integral in terms of U and T: 

^S^l±^=' .^ = e^^e--, (A.11) 
am[z^ — z^) A 

with 

e^-^ = ^T^-i ^ (A.12) 

and thus 

- (p -^z"'(P(z;x) = ^e^"'^e-''^^ , if U^O. (A.13) 
i J aA A 

The second integral in ( lA.SI) receives contributions both from the pole of A (which is at z = 2;+^ 
for U > and at z = zZ^ for t/ > 0) and from the pole of z~'^ at 2; = 0. The residue at z = z'^ 
is given by, for f/ > 

27rz— — ^ — -x) = 27r^ -i x) ^+ 

am(2;+ — 2;_ ) am(2_ — z+) 

= 2vr.^^^^±^ = 2vr. ie-^e^-^ , (A.14) 
where we have used that 

0(2:;^-x) = -0(2:+;x)- (A.15) 

Analogously one finds the residue for f/ < 0: 

^^^ cos^ + cos^- .™ = |e-^e--^ . (A. 16) 
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The residue at 2; = is given by 

1 d""-^ r</)(z;-x) 



^"^ ~ (m- 

Note that 



K 



'm 



(m 



- a A J 2=0 
a A 1 z=o 



(A. 17) 



(A.18) 



We couldn't find a compact expression for R^, nor re- write it in terms of U and T. 
Putting things together, we find 

^Im = I cosh(m[/) e-''^^ + ^ . (A. 19) 

Note that, contrary to the Page propagator, there is no cusp at U = 0. 

Some limits 

Let us first look at the limit in which — i- 0. This is achieved by taking C,' = e/2 and 6' = e, 
and letting e — )■ 0, so that ^ ae^. At r'|_ = the ip' fiber degenerates: functions of ip' are not 
regular unless they appear in the combination ^JV^e^"^' 1'^ ~ ee^'^'/^. Thus regularity requires 
that in this limit vanishes as 

Irr. ~ e'"^! . (A.20) 
One can see that the two addends in (IA.19P separately diverge as 

i cosh(mf/) e-^'"'^ ~ e-l™l , ^ ~ e"!™! . (A.21) 

However one can check with Mathematica (we could do this only up to m = 2) that the diver- 
gences cancel and that Im indeed vanishes as e'"^'. 

Suppose that one of the points, say x', is on the = surface: this happens if = r'_ or 
^' = |. One then finds that goes to a finite value. 

Finally, consider the locus \J = 0, or r_|_ — r_ = —ir'^ — r'_). This implies 6' = tt — 6. Even 
in this limit one can check that is regular (and has no cusp). 
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